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Abstract 

In the post-LHC8 world- where a Standard Model-like Higgs boson has been established 
but there is no sign of supersymmetry (SUSY)- the detailed profiling of the Higgs boson 
properties has emerged as an important road towards discovery of new physics. We 
present calculations of the expected deviations in Higgs boson couplings Kr,b, nt, kw,z, 
Kg and K^ versus the naturalness measure Aew- Low values of Aew ~ 10 — 30 give rise 
to a natural Little Hierarchy characterized by light higgsinos with a mass of /r ~ mz 
while top squarks are highly mixed but lie in the several TeV range. For such models 
with radiatively-driven naturalness, one expects the Higgs boson h to look very SM-like 
although deviations can occur. The more promising road to SUSY discovery requires 
direct higgsino pair production at a high energy e^e“ collider operating with the center- 
of-mass energy ^/s > 2/i ~ \/2 Aew^z- 
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1 Introduction 


The recent discovery of a Standard Model (SM)-like Higgs boson with mass mh = 125.09 ± 
0.24 GeV (ATLAS/CMS combined values) is a triumph of contemporary physics in that it 
provides the hrst hard evidence for the existence of fundamental scalar helds. Theoretically, 
such spinless helds are hard to comprehend due to unstable quadratic quantum corrections to 
their mass value |^. If nature is supersymmetric (SUSY), then the unwanted divergences are 
cancelled to all orders in perturbation theory thus allowing for a naturally occurring Higgs 
boson [^. Yet, so far, no sign of softly-broken weak scale SUSY has appeared at LHC |^[^. 
The growing mass gap between the weak scale- as typihed by the W, Z and Higgs boson masses 
~ 100 GeV- and the sparticle mass scale m(sparticle) >1 — 2 TeV has led to the re-emergence 
of the naturalness question: this time involving log rather than quadratic divergences. 

Recent evaluations of supersymmetric models with radiatively-driven naturalness (RNS for 
radiatively-driven natural SUSY [^) hud that for a value of Aew < 10 (30), then the gluino 
mass is bounded from above by rrig < 2.5 TeV (5 TeV). These Aew values correspond to 
'^Ew ~ (3.3%) hne-tuning respectively. In contrast, the 5a reach of LHC13 (LHC with 

= 13 TeV) for gluino pair production is estimated to be rUg ~ 1.6 TeV (100 fb“^) and 
1.9 TeV (1000 fb“^) [9,10 . In RNS models with gaugino mass unihcation, this reach can be 
extended via the same-sign diboson signature arising from wino pair production to equivalent 
2.4 TeV for 1000 fb“^ TolITTl. The upshot is that LHC13 may or may not have 


values of Trig rNj 

sufficient energy/luminosity to fully probe the entire parameter space of natural SUSYj^ 

Without a guaranteed path towards SUSY discovery, other alternatives have been explored. 
Many recent investigations promote Higgs boson prohling as a probe for physics beyond the 
SM. Since the Higgs boson h is now discovered, the goal is to measure every possible property 
of h to see if they maintain consistency with the SM or produce deviations which might point 
to new physics. These quantities include: the mass and width of the Higgs boson, its spin 
(which is essentially already determined to be spin-0 [^) and its coupling strengths to various 
SM and non-SM decay modes. The coupling strengths Ki are usually parametrized in terms 
of the SM values. Thus, Kb = ghbb/dhbbi^^)^ = 5 '/iZz/fi'/iZz(SM), etc.. Evidence for beyond 
the SM (BSM) physics would then occur from the measurement of one or more values 
{i = T, b, t, Z, W,g,'y ) to signihcantly differ (by several error bars) from the SM value of 1. 

The capability of various accelerator options to measure the Ki values has been analyzed 
15,16 and tabulated in Ref’s. 17,18 . In fact, early data from LHG8 seemed to indicate an 


enhancement in the Higgs to diphoton coupling which could have been construed as requiring 
new TeV-scale charged particles that circulate in the hyy loop 19 -21 . However, the current 


prohle of the Higgs boson is consistent with SM predictions: he., at present the Higgs boson 
appears to be “the” SM Higgs boson as no credible deviations from the SM have been found. As 
more data accrue from various collider options, the error bars on the various Higgs observables 


will tighten, and may reveal physics beyond the SM 22 24 


A particularly interesting scenario which merits investigation is that of natural supersym- 


^Prospects for LHC13 indirect searches for RNS via initial state radiation off of higgsino pair production 
reactions (monojet signal) seem pessimistic [^. Allowing for monojet radiation off of Z\Zi production, Z^ —t 
Z\l'^l~ may allow probes of the higgsino mass /r up to ~ 200 GeV assuming ~ 1000 fb“^ of integrated 


luminosity 13 . 
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metry. In fact, several previous works have already investigated this case: Ref’s. 25-28 . These 


papers all investigated models with light third generation squarks which are a consequence of 

0x2 

minimizing the “large log” contribution to the Higgs boson mass: + 

A^) In (A/msusv), where ft is the top Yukawa coupling, A is the cutoff scale, 
is the SUSY-breaking scale, and rriq^, denote the soft masses and the A-term for 

stops, respectively. The validity of this measure has been challenged in Ref’s. in that it 

sets to zero additional dependent contributions which lead to large cancellations. Alternatively, 
it is argued that the correct measure is Aew which instead requires 1. that the SUSY p term 
is comparable to the weak scale (m„eak ~ 100 GeV as typihed by the W, Z and h masses), 2. 
that is driven radiatively to negative values of magnitude comparable to ruweak and 3. that 
radiative corrections to the weak scale effective potentiaj^ (which determines the electroweak 
vacuum expectation values (VEVs) and hence the Z-boson mass mz) are comparable or less 
than m„eak- This latter condition is met for highly mixed but TeV-scale top-squarks, he., much 
heavier than values expected from large log minimization. Meanwhile, the hrst of these con¬ 
ditions implies a spectrum of four higgsino-like states Zi^ 2 , with the Z\ as a higgsino-like 
lightest SUSY particle (LSP) and candidate for dark matter. If naturalness is required as well 
in the QCD sector, then the axion solution to the strong CP problem is invoked 31 . The SUSY 
DFSZ axion model 32 not only tames the strong CP problem but also provides an elegant 


solution to the SUSY /i problem. In this class of models, the apparent Little Hierarchy as typi¬ 
hed by /i ^ m(sparticle) can be naturally generated via a radiative breakdown of Peccei-Quinn 


symmetry 33,34 


While we agree with the assessment of Ref. 35 that unnatural SUSY is likely to be wrong 


SUSY, we would disagree with the assessment that the minimal supersymmetric Standard 
Model (MSSM) is at present hne-tuned over all of parameter space. While many SUSY mod¬ 
els are indeed hne-tuned under Aew 30 , the class of SUSY models with radiatively-driven 


naturalness (RNS) remain highly natural. The reason is that the current experimental limits 
on the SUSY p parameter arise from negative searches for chargino pair production at LEP2: 

> 103.5 GeV. Roughly speaking then, /i is also > 100 GeV. This value is quite close to 
the value of mz so that we can interpret the fact that my/,mz and are all clustered near 
100 GeV as a consequence of /i^, mjj^iweak) and all being comparable to- or lighter than- 
(100 GeV)^. Meanwhile, the other soft SUSY breaking parameters can be much heavier, as is 
indicated by LHC sparticle search limits and radiative corrections to mh- 

As mentioned above, the class of RNS models predict light higgsino-like states around 
the electroweak scale. In addition, there may also be electroweak gauginos and/or heavy Higgs 
bosons below the TeV scale, which are currently less constrained as they are un-colored particles. 
The presence of such particles can in principle modify the Higgs couplings. For instance, the 
chargino loop contribution can alter the hyy coupling gh^y-y if the chargino state has a sizable 
wino component. Also, if heavy Higgs bosons have relatively small masses, the lightest Higgs- 
boson couplings deviate from those in the case of decoupling limit, i.e., the SM Higgs ones. 
The precise measurements of the Higgs couplings, therefore, may provide a way of probing the 
RNS scenario indirectly. Since forthcoming collider experiments can offer signihcantly improved 
sensitivities, it is quite important to investigate whether these experiments can actually observe 


and are given in the Appendix of Ref. [^. 
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any deviations in Higgs couplings in the case of RNS models. 

In this paper, we calculate the deviations to the Higgs boson couplings HT,b,t,w,z,g,'y in su¬ 
persymmetric models with low Aew- After a brief review of our naturalness considerations 
in Sec. in Sec. we discuss the Higgs couplings in the MSSM and in Sec. we discuss 
constraints on natural SUSY parameter space. We present in Sec. our main results of the 
values of the Ki versus electroweak naturalness measure Aew from a scan over parameters of the 
two-extra parameter non-universal Higgs (NUHM2) supergravity (SUGRA) model 36 which 
allows solutions with Aew as low as 5-10. We compare these expectations against the values 
which are expected to be probed at present and future LHC runs, and with expectations from 
the International Linear e+e“ Collider (ILC). In SUSY models with low Aew (highly natural 
models), the bulk of points give tiny deviations from SM expectation. In Sec. we show that 
the value of k,^ can be enhanced to yield deviations as high as only 2% in models with gaugino 
mass non-universality and a light wino 37 . In Sec. we compare our results against direct 
sparticle search prospects for LHC and ILC. We stress there that LHC13 has only a limited 
reach for natural SUSY. In addition, if the ILC is built initially as a Higgs factory, we ultimately 
expect from natural SUSY that the Higgs prohle will look very SM-like: any major deviation 
from the SM would likely come from a rather light spectrum of heavy Higgs bosons which 
are already highly constrained by LHC searches. We hnd that the best prospect for probing 
natural supersymmetric models remains as the direct production of higgsino pairs at ILC. In 


such a case, ILC would function as a higgsino factory and as a discovery machine for SUSY 38 


2 A natural SUSY spectrum 


Any quantitative discussion of naturalness requires the use of some measure, and several have 
appeared in the literature. Before proceeding, however, we note the observation that some 
measures can be mis-applied by claiming large opposite-sign contributions to observables of 
dependent quantities: these mis-applications lead to over-estimates (29| of electroweak hne- 
tuning]^ 

When 


To avoid such pitfalls, any naturalness measure should obey the hne-tuning rule 30 


evaluating fine-tuning, it is not permissible to claim fine-tuning of dependent guantities one 
against another. 


2.1 Electroweak scale naturalness 

The relationship between the weak scale m„eak and SUSY parameters arises from minimizing 
the MSSM scalar potential. One is led to the relation 

ml _ + ^i- + K) fi 2 ..n 

2 tan^ /9 — 1 ^ ’ 

for the Z mass mz, where and denote the 1-loop corrections (expressions can be found in 
the Appendix of Ref. j^) to the scalar potential, and are the Higgs soft masses, and 

^For example, if an observable O is expressed as O — O -\-h—b where b is large, then O may look fine-tuned. 
In this trivial example, combining dependent contributions then cancels the would-be source of fine-tuning. 
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tan/3 = {Hy)/{Hfi) is the ratio of the Higgs VEVs. SUSY models requiring large cancellations 
between the various terms on the right-hand-side of Eq. Q to reproduce the measured value of 
m| are regarded as unnatural, or hne-tuned. In contrast, SUSY models which generate terms 
on the RHS of Eq. Q which are all less than or comparable to mweak are regarded as natural. 
Thus, the electroweak naturalness measure Aew is dehned as |8|[40] 

Aew = niax|each additive term on RHS of Eq. 0 |. (2) 


Including the various radiative corrections, over 40 terms contribute, 
corrections, and taking moderate-to-large tan/3 > 3, then n7.|/2 ~ 
criterion for naturalness is that at the weak scale 


Neglecting radiative 
— so the main 


-mt 


and /i 


^ ~ mi 


41 


The value of 


m'jj (where niA ~ miif^(weak) with rriA being the mass of the CP-odd Higgs boson) can lie in 
the TeV range since it is suppressed by 1/tan^/3. The largest radiative corrections come from 
the top squark sector. Requiring highly mixed TeV-scale top squarks minimizes 2 ) whilst 
lifting the Higgs mass rrih to ~ 125 GeV [^. 

Some virtues of Aew include 1. that it is model independent so that any model generating 
the same weak scale spectrum will have the same naturalness value and 2. it obeys the hne- 
tuning rule. It is also predictive: since |/r| ~ mz, it implies a spectrum of four higgsino states 
Z 12 and all lying not-too-far from mz- /i ~ 100 — 200 GeV. While many models are indeed 
highly hne-tuned under Aew 30,42 , the NUHM2 model and its generalizations admit Aew 
values as low as 5 — 10 leading to just 10 — 20% electroweak hne-tuning. The models with low 
Aew ^ 30 are regarded as natural. 


2.2 Higgs mass fine-tuning 


An alternative measure comes from Higgs mass hne-tuning 43 45 . The light Higgs mass for 
tan /3 > 3 is given by 

ml ~ -2{m^^ (A) + 5m]j^ + , (3) 

where the largest contribution to includes divergent logarithms of the ehective theory cut- 
oh scale A which is commonly taken to be A ~ tugut ~ 2 x 10^® GeV in gravity-mediation. By 
neglecting gauge terms and setting mk\j^ to zero, a one step integration of the renormalization 

group equation (RGE) for m\j^ leads to ~ — ^(mg^+m^g-l-Af) In (A/rngusv)- Requiring 

< ml then implies the existence of three third generation squarks ti^ 2 , with mass less 

It also leads to claims that SUSY is hne-tuned at the per-mille 


5m?jT 

than about 600 GeV 45 


level 46 . 


The problem with this measure is that m 


and Smir 

J^u 


are 


not independent!^ In fact, the 


larger the value of (A) becomes, then the larger becomes the cancelling correction 47 
combining the dependent terms (m^ (A) -|- 5m‘jj ), then instead one is lead to requiring unau 
and the weak scale value of are both comparable to ml- Even if m\j^{K) lies in the 
multi-TeV range, it can be driven to small negative squared values at the weak scale m^eak via 
the same radiative corrections that drive electroweak symmetry breaking in SUGRA models. 


^This is different from the case of the SM Higgs mass fine-tuning where the tree level mass and quadratic 
divergences are independent. 
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After combining dependent contributions to then a low Higgs mass fine-tuning implies the 
same general consequences as those of low Aew- 


2.3 BG fine-tuning 

The BG measure |48|-[50] is dehned as 


Abg = maxi 


d log 
d\ogpi 


(4) 


where Pi are fundamental (usually high scale) parameters of the model labeled by index i. To 
evaluate, we start with the weak scale relation 


m| ~ —2p^(weak) — 2m^^(weak) , 


(5) 


and express the right-hand-side in terms of fundamental high scale parameters. A pitfall can 
occur in what constitutes high scale independent parameters. Since p hardly evolves during the 
renormalization group (RG) flow- it only receives the wave-function renormalization thanks to 
the non-renormalization theorem of the superpotential- then p(weak) ~ p{A). On the other 
hand, evolves greatly: indeed, it must be driven through zero to negative values by 

the large top quark Yukawa coupling ft for electroweak symmetry to be broken radiatively. 
Semi-analytic solutions to the RGE allows m|^^(weak) to be evaluated as a large sum of 
contributions from various high scale soft parameters: some positive and some negative. In the 
case of gravity-mediation, however, for any particular hidden sector the high-scale soft terms 
are all calculable as multiples of the gravitino mass 7713/2 [^. If we vary 7773 / 2 , the soft terms 
all vary accordingly: i.e. they are not independent in SUGRA models. By combining the 
dependent soft SUSY breaking terms, then the Z mass can be expressed as 

ml ~ -2/i^(A) - 0777^/2 , (6) 

with a being a certain proportionality factor dependent on each soft mass spectrum. Using 
Eq. ([^- and since p hardly evolves from A to TTr^eak^ we have — 2mjj^{weak). Even if 

7773/2 is large (as implied by LHG8 limits for gravity-mediation), then one may still generate 
natural models if the coefficient a is small. Under the combination of dependent soft SUSY 
breaking terms, then low Abg implies the same as low Aew^ that p ~ 777„eak and that is 
driven to small and not large negative values. 


3 Higgs couplings in the MSSM 

In this Section, we briefly review the Higgs couplings to SM particles. In the MSSM, the 
lighter of the two GP-even Higgs mass eigenstates is typically a SM-like Higgs boson but with 
properties differing from the SM case depending on the mixing angle a. In general, the Higgs 
boson couplings to vector bosons {W and Z bosons) are simply determined by a and (3 while 
couplings to fermions have contributions from loop corrections as well. On the other hand, the 
dimension-hve couplings of Higgs to diphoton and to digluon are generated at one-loop order. 
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Note that also in the SM these couplings are induced at loop level. For this reason, these 
couplings can be rather sensitive to the SUSY effects. 

In the MSSM, the SM-like Higgs boson is usually the lighter eigenstate of the mass matrix 

+/i2 + m|(l-2cos2/3)/2 + 2/^2 + m|) sin2/?/2\ 2 

^ + 2p2-I-77i|) sin2/d/2 -p -|-2 cos2/3)/2 ) ^ 

(7) 

in the basis of the weak eigenstates of the CP-even neutral Higgs helds The radiative 

corrections to the Higgs mass-squared matrix are included in It is worth noting that we 

use Higgs soft masses and /r in order to directly compare it with hne-tuning argument in Eq. (|^. 
The conventional form of mass matrix is obtained if one uses the relations, 


uIa = m 


„ + + 2/x' 


H: 

--COS 2d 

COs2/3 y 2 ^ ^ 

("'L + d' + cos 2(3 ) , 


( 8 ) 


from the electroweak symmetry breaking conditions. The mass-squared matrix is diago¬ 
nalized by the mixing matrix]^ 


( h\ _ f cos a sina!\ 
\Hj sina cosay 


The vector boson couplings are simply given by 


9hvv = glvv sin(a -F (3) for V 


w,z. 


( 10 ) 


In the decoupling limit where rriA 3> the mixing angle a follows the relation, a + (3 ~ 7r/2. 


This decoupling behavior can be clearly seen in the approximate formula 52 


/ '^7 

cos(a + (3) = —- 


( 11 ) 


Note that if the radiative corrections in Eq. 0 are sub-dominant, then cos(a -|- /3) < 0 since 
sin 4/3 < 0. This determines the direction of the deviation in the Higgs-fermion couplings, as 
we will discuss below. The above equation reads 


sin(a -1- /3) ~ 1 — 


1 

2 


cos2(a -1- (3) 


1-0 



( 12 ) 


From this relation, one can easily see that cos(a -|- /3) —)■ 0 and sin(a -f /3) —)■ 1 as —>■ 00 . In 
addition, we note that the deviations in the gauge boson couplings from the SM ones rapidly 
decrease in the large tua limit since they are suppressed by a factor of m%(m\. Thus, we expect 
that the vector boson couplings are almost SM-like, as will be actually seen in Sec. 


®In the notation of Ref. 


39 


the Higgs mixing angle a is the negative of other conventions. 


6 










For the case of fermion couplings, we need to consider an effective Lagrangian|^ for the 
Yukawa couplings. Including SUSY threshold loop-corrections, the low-energy effective Yukawa 
terms below the SUSY breaking scale are written as |54|[55] 

- -Ceff = ih + 6fb)bRHdQ3 + AhbRH:Q, + {ft + 6ft)fRH^Q3 + AftiRH^Q^ + h.c., (13) 

where 6fb and 6ft represent the radiative corrections to the tree-level bottom and top Yukawa 
couplings fb and ft in the MSSM superpotential, respectively, while Afb and Aft are loop- 
induced non-holomorphic Yukawa couplings. Notice that these (non-logarithmic) radiative 
corrections are generated by the SUSY breaking effects; in the SUSY limit, the vertex cor¬ 
rections to the Yukawa couplings vanish because of the non-renormalization theorem. These 
radiative corrections modify the relations between the fermion masses and the corresponding 
Yukawa couplings as 


mb = fbV cos (3 ^ tan (3^ = fbV cos (3 {1 + Ab) , 

mt = ftv sm 13 (^1 + ^ cot 13 '^ = ftv sin 13{1 + At) , 


where u ~ 174 GeV denotes the Higgs VEV and are given by 

m? mA -n -1 

bi 


Ab ~ 

At ~ 




2 


tan l3 




(14) 

(16) 


(16) 

(17) 


Here, a* denotes the strong gauge coupling constant. M 3 and are the masses of gluino 
and sbottoms, respectively. The loop function is dehned by 


I{a,b,c) = 


abln{a/b) + bcln{b/c) + caln{c/a) 
{a — b){b — c){a — c) 


(18) 


This function is of the order of l/max{a^, 6 ^, c^} [5^. Notice that for large tan/?, Ab — 
Afbtan(3/fb since the non-holomorphic correction is enhanced by tan/?. The non-holomorphic 
correction to the top Yukawa coupling is, on the other hand, suppressed by tan /?, and thus At ~ 
bft/ft in this case. The modihcations in the relations (14) and (15), as well as the deviation 
from the decoupling limit characterized by Eq. 0 . change the fermion-Higgs couplings from 

®In our analysis, the low-e nerg y effective theory is matched onto the full MSSM at the scale of msuSY = 


in the DR-scheme 53 , with ^ being the masses of stops. 
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the SM ones. We hav4l] 


Qhhh 

= qH 

Qhtt 

= 9m 


. / cos(a +/3) f o A ^ ^ 

sm(Q! + /9)---< tan (3 — Ai, cot [3 + (tan [3 + cot (3) —— 

1 + Afe ( fb 


sin(Q; + (3) + 


cos(a + /?) , A A ,, , a o\^ft 


1 +At 


(1 + At) cot /3 — (1 + cot^ (3) - 


ft 


, ( 21 ) 
( 22 ) 


Here we express the latter equation in terms of At and Aft since At ~ 3ft/ft- One can also 
obtain a similar relation for the Higgs coupling to tan lepton, g^rr^ by replacing Ab with A,- 


and Sfb/fb with Sf^/fr in Eq. (21). In this case A,- is dominantly induced by the wino-higgsino 
loop diagram 1 ^: 

* ’ 3q( 

Ar ^ —^M 2 /xtan/?/(m~ ,M 2 ,/i^) , (23) 

OTT 

where a 2 is the SU(2)i gauge coupling constant, M 2 is the wino mass, and is the left-handed 
third generation slepton mass. Here, we have dropped the sub-dominant bino contribution for 
brevity. From the above equations, it is found that the deviation from the SM couplings is 
proportional to m|/m^ and therefore becomes quite small in the large limit. In addition, 
we see that the deviation in the bottom and tan couplings is enhanced by tan f3, while that in 
the top coupling is not. Thus, the bottom and tan couplings are more appropriate to probe 
the SUSY effects than the top coupling, as will be shown below. Moreover, we note that as 
long as the radiative corrections are moderate, the bottom/tau (top) coupling is always larger 
(smaller) than the SM one as cos(a - 1 - /d) < 0. This feature is also found in the analysis given 
in Sec. [5l 

For the various Higgs couplings to gg pairs, we use the standard expressions including 
quark and squark triangle diagrams as given in Ref’s [^. We hnd that the SUSY effect 
on the effective gluon coupling is dominantly given by the stop contribution, which can be 


approximately expressed as 58 




A mt 
rs-' 1 _L -L 

9 - ^ + 4 


+ 


m 


ti 


1 {At — p cot /d)^ 
mi mi mi 

C2 Cl C2 


(24) 


This expression shows that Kg <1 occurs only if the stop mixing is sizable. As discussed above, 
the RNS scenario with mh — 125 GeV requires large stop mixing, and thus we expect that the 
gluon coupling can be smaller than the SM one in this scenario. 

For the Higgs couplings to 77 , we use standard expressions including quark, lepton, squark. 


slepton, and IV ^2 loops 57 . In the RNS models, higgsinos lie around the electroweak 

scale, and thus may give rise to a considerable contribution to the 77 coupling if the higgsinos 
well mix with winos to have a sizable coupling with the Higgs boson. This can be achieved 
when wino has a relatively small mass. We discuss this possibility in Sec. 


^Here, we have used the identities 


sin a 

cos P 
cos a 

sin/3 


sin(Q! + P) — tan /3 cos(a + /3), 


sin(Q; + P) -\- cot/3cos(a -I- /3). 
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(19) 

( 20 ) 




















4 Constraints from low energy observables 

We next discuss some low energy observables that constrain the low niA region in our scanned 
results as shown in the next Section. As discussed in the previous Section, Higgs couplings to 
vector bosons and fermions suffer deviations mainly from Higgs mixing so that large deviations 
in the Ki are mainly expected when is small. In such cases, however, loop-mediated con¬ 
tributions to B decays and also tree-level contribution mediated by charged Higgs also become 
larger. In this case, H-decay observables can constrain this portion of parameter space. 


4.1 BR{B X,-f) 


For our evaluation, we use the NLO SUSY calculation from Ref. 59 as included in Isatools. In 


the MSSM, the two major SUSY contributions come from chargino-stop loops and also from 
the charged Higgs-top loop. In the large tan /3 regime, these are approximately given by 

BR(B ^ X,A\n-,i ■ ■ " 


oc At tan (3- 


BR{B —)■ XsX\H±,t cc 


n(l -f Ab) 
mbift cos/3 - A ft sin/3) 
V cos (3{1 + Ab) 




(25) 

( 26 ) 


where / and g are loop functions 60 . In order to lift the Higgs mass to 125 GeV, it is 


normally required that stop masses are of order TeV although one of them can be below TeV 
if the maximally-mixed stop scenario is considered. In such a case, then the chargino-stop loop 
contributions are usually small. For the small case- where Higgs mixing can sizably affect 
vector boson and fermion couplings- the light charged Higgs ~ m\ + ruw) can make a 

sizable contribution to the decay width. In this case, the small rriA region is typically excluded 
by the BR{B —)■ A^y) measurement. 


4.2 BR{Bs -A ) 


The Bs —)■ gAgT decay is induced by flavor changing interactions of Higgs bosons, h, H, and 
A. The flavor changing couplings of Higgs with b- and s-quarks are generated by similar loop 
processes as those for Ab and At. The physical discussion and calculation details are provided 
in Ref. 61 . Since ghbs cc cos(a + (3)/ cos^ (3, guts oc — sin(Q! + {3)/ cos^ (3 and gAbs oc —1/ cos^ /?, 


the dominant contributions are from H and A mediated processes. The overall branching ratio 
is then given by 
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BR{Bs —)■ g'^g ) cx taxr (3 — -r- 


m, 


(27) 


Hence small nriA and large tan (3 enhance BR{Bs 
stringently restrict such a parameter region. 


-A g^g ), and thus the experimental bounds 


5 Results for RNS in NUHM2 

In this Section, we explore the Ki values which are expected in SUGRA GUT models with 
low hue-tuning Aew- We will adopt the two-extra parameter non-universal Higgs model 
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(NUHM2) as a template. The parameters in this model are given by 


mo, mi/ 2 , Ao, tan/3, p, m^. 


( 28 ) 


In the above set, tuq is the GUT scale value of the common soft mass parameter for matter 
scalars, mi /2 is the unihed gaugino mass, Aq is the unified trilinear soft term, tan/3 is the usual 
ratio of VEVs and fi and niA are the weak scale values of the superpotential fi parameter and 
the pseudoscalar Higgs mass. These latter two determine the weak scale values of m|^^ and 
m|^^ from the scalar potential minimization conditions. The GUT scale values of m|^^ and 
are determined by RG evolution and are then in general non-universal with the matter scalar 
mass mo. 


over the 

parameter space 

mo 

: 0-20 TeV, 

mii2 

: 0.5-2 TeV, 

-3 

A 

O 

o 

A 

h 

: 0.1 - 1.5 TeV, 

rriA 

: 0.15 - 20 TeV, 

tan (3 

: 3-60, 


(29) 


and generate sparticle mass spectra using Isajet 7.84 62 which contains the Isasugra subpro¬ 
gram. The range of /i covers only positive values; the physical results are very similar in the 
case fi < 0 except that ABr(6 —)■ sy) limits are more constraining for low m^ and hence we get 
smaller deviations in the Ki values. The major difference with negative fi is which will be 
discussed in Sec. El 

We require of our solutions that: 


• electroweak symmetry be radiatively broken (REWSB), 

• the neutralino Zi is the lightest MSSM particle. 


• LHG search bounds on rrig and rriq in mSUGRA are respected [^[^, 

• -2.3 X 10-9 < ABr(R, ^ p+p-) < 0.6 x 

• —3.6 X 10“^ < ABr(6 —)■ sy) < 9.2 x 10“® 

• rrih = 125 ± 2 GeV. 

Here, we have taken a ±2 GeV error range in the Higgs mass rrih to reflect the theoretical 
uncertainty of the computation]^ The lower bound on the Higgs mass rules out the region 
where mo ^ 0.4 TeV. The upper bound on the branching ratio of 6 —)■ sy decay removes the 
bulk of parameter region where m^ < 0.3 TeV. 

®We implement the RG-improved 1-loop effective potential calculation of rrih which includes leading 2-loop 
terms in Isajet. 



the light chargino mass obeys the model independent LEP2 limit, m^^ > 103.5 GeV 
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Figure 1: Kr vs. Aew from scan over NUHM2 parameter space with rrih = 125 ± 2 GeV and 
LHC Higgs and sparticle mass constraints. The current LHC reach (gray-dashed) and future 
reach of LHC (green-dashed) and ILC (red-dashed) are shown (see the text for details). 


Our hrst results are shown as Kr vs. Aew in Fig- 0 Here, the dots are color coded as to the 
value of niA, with blue indicating < 0.5 TeV, orange is 0.5 TeV < uia < 0.75 TeV, purple is 
0.75 TeV < niA < 1 TeV and black is > 1 TeV. As discussed in Sec. this coupling has a 
large deviation from unity if cos(a + /3) is sizable, which occurs when is comparable to ruh. 
Thus, the magnitude of Kr follows the mass values for the heavy Higgs eigenstates: a value of 
Kr ^ 1 when the heavy Higgs eigenstates decouple. Furthermore, in a wide range of parameter 
space Kr is larger than unity, as discussed in Sec. In the RNS model with low Aew, Ihl ~ 


and \mHj ~ ruA- Since enters Aew as ~ tan^/3, then rather large values of mA 

are consistent with low hne-tuning. Upper limits on rriA have been found in Ref. 64 where 
tiia <5 — 8 TeV for Aew < 30 (the exact upper bound depends on tan/?). Thus, the bulk of 
points with relatively large ttia and low Aew are expected to give only slight deviations from 
the SM hrr coupling. The points with large deviations occur for low niA, and in fact there are 
already tight constraints from LHC on gg ^ h, H, A ^ for SUSY in the ruA vs. tan f3 
plane |^. These LHC heavy Higgs search constraints need revision for RNS SUSY since in 
RNS the requirement of rather light higgsinos means the heavy Higgs bosons dominantly decay 
to charginos and neutralinos 64 rather than to SM modes such as r+r“. We do impose the 
ATLAS rriA vs. tan (3 constraints in the case where mA,H < 2/r. 

In the plot, we also show the current reach for Kr from LHC8 as the gray dashed line at 
Kr ~ 1.3 and the future reach of high luminosity LHC13 (HL-LHC) with 3000 fb“^ and ILC500 
in the green- and red-dashed lines, respectively 17,^. From current reach of LHC8, we can 
conclude that the LHC experiment has already disfavored mass spectra with mA < 300 GeV. 
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Figure 2: vs. Aew from scan over NUHM2 parameter space with m/j = 125 ± 2 GeV and 

LHC Higgs and sparticle mass constraints. 


Furthermore, it turns out that the high-luminosity LHC as well as the ILC can probe ruA ~ 
1 TeV, which can be clearly seen in the inset of the hgure where a magnihed view of Kr very 
close to 1 is shown. While future colliders can probe much of the parameter space with low 
Aew, a large chunk with multi-TeV values of would look very SM-like. 

In Fig. 1^ we show Kb vs. Aew- The locale of the dots is nearly the same as for the v plot 
since Kb ~ Kr- The main difference occurs in the current and future collider reach for Kb- Here, 
HL-LHC is expected to probe a 4% deviation while ILC500 can probe a 1% deviation. While 
these reach values probe a large fraction of parameter space with < 1 TeV, there are a 


number of natural models which predict quite small deviation in Kb- 

In Fig. 1^ we show the values of Kt vs. Aew- As discussed in Sec. this coupling is expected 
to suffer hardly any deviation from the SM value, since there is no tan (3 enhanced effect in this 
case. In addition, the projected experimental probes are much more limited since the h ^ tt 
decay mode is kinematically closed. The value of Kt must be extracted from hts to the Higgs 
production coupling hgg which includes a top-quark loop in the case of LHC, and also to tih 
production in the case of LHC and ILC. Here, it is expected that ILC500 may probe to the 
2.5% level {Kt ~ 0.975) once a/s > 2mt -|- ruh- 

In Fig. 1^ we show the values of kz vs. Aew- In this case, the value of kw,z is expected to be 
close to 1 since the deviation is suppressed by m^jm\ as shown in Eq. (12). On the other hand, 
HL-LHC can probe kz to ~ 2% precision via h ^ ZZ* ^ M decays and ILC can probe to 
sub-percent precision since h is dominantly produced via Higgsstrahlung: e^e“ Z* ^ Zh- 
Even so, the bulk of points with low Aew have only tiny deviations from 1 and so in this 
channel one expects the h to look highly SM-like for RNS SUSY. 
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Figure 3: vs. Aew from scan over NUHM2 parameter space with ruh = 125 ± 2 GeV and 

LHC Higgs and sparticle mass constraints. 


As we have discussed so far, the deviations in the fermion and gauge boson couplings are 
mainly due to the effects of a sizable cos(a + (3), which occurs if is relatively light. Such 
effects are, however, also induced in the two-Higgs doublet models. To conhrm the presence of 
SUSY effects, therefore, it is desirable to see the contribution given by other particles than the 
Higgs bosons. To that end, we consider the following quantity discussed in Ref. 


52 


Kb 


Kt 


Kt - Kb 


^ Afe , 


(30) 


where we have kept only the tan/^-enhanced terms and used the fact that |At-| is rather small 
since it is induced by the electroweak gaugino loop diagrams. In Fig. we plot this quantity 
vs. Aew with color coding in accord with mi^. It is found that a sizable value of A;, is 
expected in most of parameter points. Therefore, we may extract even the information of the 
sfermion/gaugino sector via the precise measurements of the fermion couplings. However, we 
also note that the value of the quantity is found to be relatively small in the small Aew region. 
This challenges the extraction of Af, in the natural SUSY scenario. 

In Fig. we show the value of Kg vs. Aew- In this plot, the dots are color coded as 
to the value of with blue indicating < 1.5 TeV, yellow is 1.5 TeV < mq < 3 TeV, 
green is 3 TeV< mq < 4 TeV and red is mq > 4 TeV. Here, the hgg coupling proceeds 
from triangle diagrams including quarks for the SM case and also squarks for the SUSY case. 
Thus, we expect large deviations from the SM coupling if squarks are far lighter than the TeV 
range. Since we require rrih ~ 125 GeV, then we implicitly require TeV-scale highly mixed top 
squarks which provide a sufficiently large radiative correction to rrih- Usually the top squarks 
are amongst the lightest squarks since their masses are suppressed by large top-quark Yukawa 
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Figure 4: kz vs. Aew from scan over NUHM2 parameter space with mh = 125 ± 2 GeV and 
LHC Higgs and sparticle mass constraints. 



• rrifi < 0.50 TeV 

• 0.50 TeV <m^< 0.75 TeV 
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effects in RG running, and also by large mixing. Furthermore, the NUHM2 model should obey 
well the LHC 8 constraints on the minimal SUGRA model (mSUGRA) so that rriq >1.8 
TeV. Thus, we do not expect squarks well below the TeV-scale and therefore large deviations 
in the Kg coupling. While some points with low have deviations of several %, which can 
be probed by HL-LHC via the overall s-channel Higgs production rate (j{gg —)■ h), the bulk of 
points with a decoupled in the TeV-range tend to have deviations of less than a percent. 
These deviations will be hard to access by either HL-LHC or by ILC. Note that most of the 
parameter points predict Kg < 1 . As discussed above, this can happen only when there exists 
a large left-right mixing in the stop mass matrix. Since this large mixing is a typical feature of 
the RNS models, the reduction in the hgg coupling can be regarded as a distinctive prediction 
in the RNS scenario. 

In Fig. we show the value of k^ vs. Aew- Color coding is the same as in Fig. In the 
SM, the K^ coupling proceeds via triangle diagrams involving charged particles which couple to 
the Higgs: the gs, and W^. Among them, top quark and W boson give rise to the dominant 
contributions. In the case of SUSY, then there are additional loops containing squarks, sleptons, 
charginos and charged Higgs bosons. As in the case of Kg, large deviations are obtained in the 
light stop region which also coincides with the small p region (with light charginos). Moreover, 
for mH± ~ mA small (large Higgs mixing) then the hH^H~ coupling can lead to deviations in 
K^ as in two Higgs doublet models. For light charginos, then the hW\Wi coupling can be large 
and also contribute. Thisj»uphng is proportional to higgsino-times-gaugino components of Wi 
and so in the case where Wi is nearly pure higgsino or wino, the coupling is smaller. From the 
plot, we expect deviations in k^ < 1%. 

Even though the h —)■ 77 branching fraction is small, the LHG gg ^ h production cross 
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Figure 5: {Kb—Kr) / {Kt—Kb) vs. Aew from scan over NUHM2 parameter space with mh = 125±2 
GeV and LHC Higgs and sparticle mass constraints. 
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0.05 


0.00 
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section is large and the 77 signature is robust. For comparison, the reach in K<y of HL-LHC is 
shown which extends to the 2% level. This is not enough to access the bulk of low hne-tuned 
RNS models. The small h —)■ 77 branching fraction limits the ILC capability to probe this 
loop-induced coupling. ILC500 is projected to probe values of to the 8 % level. 0 


6 Natural SUSY with light wino 

The results from the previous Section were evaluated in the NUHM2 model which assumes 
gaugino mass unihcation: Mi = M 2 = M 3 at the GUT scale so that M 3 ~ 7Mi and M 2 ~ 2Mi 
at the weak scale due to RG evolution. Then the LHC limit (that nig > 1.3 TeV from the 
mSUGRA cascade decay analysis) translates roughly to M 2 > 350 GeV and Mi > 175 GeV. 
This means for RNS SUSY with low /x ~ 100 — 200 GeV that the light Wi which circulates in 
the h 77 loop is mainly higgsino-like and has somewhat suppressed couplings. The hyy coupling 
can be increased in models with non-universal gaugino masses where nig can remain above the 
LHC 8 bound, but now M 2 and Mi can be much lower resulting in natural SUSY with either a 
wino-like or bino-like LSP 


37 


In the RNS case with non-universal gaugino masses and a lower value of M 2 , then the Wi 
can be a wino-higgsino admixture. Such a mixed chargino enhances its coupling to the Higgs 
bosonj^ h which depends on a product of gaugino times higgsino components. 

To show a case with maximal in RNS, we plot in Fig. [^the value of vs. M 2 along an 
RNS model-line with parameters mo = 5 TeV, mi /2 = 0.7 TeV, Aq = —8 TeV, fi = 200 GeV 


^See p. 178 of Ref. 


39 


15 












1.00 


0.99 


O) 







0.98 


0.97 


0.96. 


* m- <1.5 TeV 

, 1.5 TeV < m- <3 TeV 
, 3 TeV < m- <4 TeV 

* m- >4 TeV 

h 


0 


20 


40 


60 


— ILC500 

— HL-LHC 


80 


100 


^EW 


Figure 6: Kg vs. Aew from scan over NUHM2 parameter space with rrih = 125 ± 2 GeV and 
LHC Higgs and sparticle mass constraints. 


and niA = 1 TeV. We abandon gaugino mass unification and instead allow M 2 to vary. We 
also plot remits for several tan /? values. For the case shown, as M 2 decreases fromjts universal 
value, the PFi becomes more of a mixed wino-higgsino state and the coupling hWiWi increases. 
Correspondingly, k^ increases. The maximal k-^ reaches ~ 1.03 for lower values of tan f3 and 
for M 2 ~ 150 GeV with Aew ~ 10. Such a large value of k^ should be accessible to HL-LHC 
as shown by the green dashed line. Larger values of tan /3 > 40 (Aew ~ 50) are excluded by 
Bs —)■ iJi~ constraint since Br(Hs —)■ oc tan®/3 as stated in Eq. (27). 

It is also interesting to see that negative p makes K-y smaller than unity. If stops are as 
heavy as a few TeV, which is required to obtain the 125 GeV Higgs mass, main contributions 
to Higgs-to-diphoton decay come from chargino loops, so < 1 means that chargino loop 
contributions destructively^ interfere the dominant W boson loop contribution. It is simply 
understood from the hWW coupling, which is given by^°' 


9hw,m - 92 sign(/i) 
9hW2W2 - 92 sign(/i) 


m^{M 2 cos /d -|- /isin (5) 


Mi 




m^(M 2 sm(3 + cos /3) 


M|-p2 


5 


(31) 

(32) 


where m\ S> m\, |M| — /i^| 3> and \M 2 /fi\ < tan/?. Here we assume that hFi is mostly 
higgsino-like and kF 2 is mostly wino-like. The chargino-Higgs couplings flip their sign when the 
sign of /i is flipped, and thus chargino contributions can be either constructive and destructive 


^°See Sec. 8.3 and 8.4 of Ref. 
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for complete formulae. 
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LHC Higgs and sparticle mass constraints. 
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depending on the sign of p. In order to avoid chargino LSP, for /r < 0 we set Mi = 100 GeV. 
K.y can show about 2% dehcit for small tan f3 and M 2 , and it approaches to the SM value as M 2 
increases (black curve). 

If a deviation in is actually observed at the LHC, this may indicate the presence of 
a light chargino with sizable coupling to the Higgs boson. Such a light chargino should be 
within the reach of direct production at the ILC. Moreover, a large coupling to the Higgs 
boson implies a large neutralino-nucleon scattering cross section. Although we expect that 
only a small portion of dark matter energy density is occupied by the Zi LSP since such a light 
higgsino-like neutralino in general results in a small relic abundance, it is found that future 
dark matter direct detection experiments can probe the Zi LSP in this case, and thus provide 
a way of examining RNS models 
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7 Conclusions: 

In this paper we have presented expectations for possible deviations in Higgs couplings that are 
expected in SUSY models with radiatively-driven naturalness. Such models with low Aew ^ 30 
are natural in the electroweak sector and, if augmented with a Peccei-Quinn sector, are natural 
in the QGD sector as well. Models with a SUSY DFSZ axion also admit an elegant solution 
to the SUSY fi problem. Such natural SUSY models are consistent with squark, slepton and 
gravitino masses in the multi-TeV range which admits a solution to the gravitino problem |67p8] 
and at least a partial decoupling solution to the SUSY flavor and GP problems (^. These 
models are rather simple extensions of the SM and may even be regarded as more conservative 
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Figure 8 : vs. M 2 along the RNS model line for various values of tan jS. The reach of HL-LHC 

is shown by the green dashed line. 


than the SM in that they contain solutions to the gauge hierarchy and strong CP problems. 
Thus, every avenue for their verihcation must be explored. Here, we examined the case of Higgs 
boson prohling. 

Our results may be summarized as follows. 


• Substantial deviations in Hj. and Kf, may be expected for RNS SUSY but mainly in the 
case where niA is rather light leading to significant mixing in the scalar Higgs sector. 
However, since can extend into the multi-TeV range at little cost to naturalness (due 
to tan^/3 suppression of the term including in Eq. ([^) these deviations may well he 
below the reach of HL-LHC and even ILC500. 


• Tiny deviations to Kt are expected. This coupling is also difficult to measure unless one 
has a linear collider with ^/s > irth + 2 mt. 

• Tiny deviations are expected in usually below the 0.5% level. 


Some deviations can occur in the Kg coupling, but mainly for anomalous cases with very 
light top squarks U. However, light top squarks generally lead to large deviations in 
BF{h —>■ 57 ) and also have recently been tightly constrained by top-squark pair produc¬ 
tion searches at LHC 8 70,71 . Except for such cases, most RNS predictions for Kg he 


below the reach of HL-LHC and ILC500. 


• Small deviations in k^ are expected- usually at the sub-0.5% level- below the reach of 
HL-LHC and ILC500. However, in models with non-universal gaugino masses where the 
light chargino becomes a wino-higgsino mixture, then k-^ may increase to the 1-3% level. 
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Figure 9: (a) Plot of Aew vs. rrig from a scan over NUHM2 parameter space, (b) Plot of Aew 
vs. from a scan over NUHM2 parameter space. We also show the curve Aew = (^1/2) 

and the reach of various ILC energy options for higgsino pair production. 


To summarize: except for unusual cases (highly mixed Higgs sector with low tua, anoma¬ 
lously light stops soon to be excluded by LHC or highly mixed charginos) natural SUSY predicts 
minimal deviations from a SM-like portrait of the light Higgs boson. Given this situation, it is 
useful to compare these indirect search methods against the direct search for natural SUSY at 
LHC and ILC. 

A direct search for gg production at LHC13 with 1000 fb“^ can reach up to rrig ~ 2 TeV 10 . 
This approximately covers Aew < 7 as seen in Fig. |^a). The LHC13 1000 fb“^ reach for gg 
production is also shown in terms of Aew by the brown histogram of Fig. LHC13 can also 
search for light higgsino pair production pp —)■ where Z 2 —)• pT Z^. Since the dimuons 

tend to be rather soft (due to the small mass gap) a trigger on hard jet radiation 

from the initial state is needed [o. The reach of various LHC13 options for p~^p~j + 
production is also shown in Fig. 

The most direct test of SUSY naturalness occurs via the direct search for higgsino pair 
production at an e+e“ collider with ^/s > 2p. Such a machine would be a higgsino factory 38 


in addition to a Higgs factory. The value of is plotted versus Aew in Fig- |^b) which 


exhibits the tight correlation where -y/s > and where 


~ /i. Since Aew ~ 

then ILC probes directly values of Aew according to ^/s ~ y/2K^mz- From the plot, we see 
that ILC500 makes a complete probe of Aew < 15 and ILCIOOO probes Aew < 55. 

we show the reach in Aew of prospective experiments. ILCIOOO can see the 


In Fig. 10 


entire RNS parameter space whereas LHC14 and TLEP can probe only a portion of it. Light 
higgsinos should ultimately be detected at ILC with -y/s > 2/i. 
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